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Abstract 

Thin films of binary mixtures that interact through isotropic forces and di- 
rectionally specific "hydrogen bonding" are considered through Monte Carlo 
simulations. We show, in good agreement with experiment, that the sin- 
gle phase of these mixtures can be stabilized or destabilized on confinement. 
These results resolve a long standing controversy, since previous theories sug- 
gest that confinement only stabilizes the single phase of fluid mixtures. 
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The phase behavior of mixtures in thin films can be significantly different from the 
bulk 0-0. Fisher et al. theoretically considered an incompressible fluid mixture which 



interacted with isotropic nearest neighbor forces and showed that confinement only stabilized 
the single phase 0-0]. This prediction is robust, and cannot be altered by the magnitude of 
the interaction between the fluids and the walls. 

Experiments, however, show that, depending on the surfaces, confinement can either 
stabilize or destabilize the single phase relative to the bulk [^-[1^. Since the model of 



Fisher corresponds to the simplest case of a binary fluid with isotropic nearest neighbor 
forces, a model with more complex interactions could rationalize the experiments |^. We 
have performed Monte Carlo simulations on mixtures which interact with isotropic nearest 
neighbor interactions and directionally specific bonds. The competition between the unfa- 
vorable dispersive interactions, and the favorable "hydrogen bonds" (HB) leads to closed 
loop phase diagrams in the bulk. While the bulk behavior of HB systems has been well ex- 
plored [|I^-[T^ , their behavior in confined geometries has not been studied theoretically. We 



find that, when these systems are confined between plates which only interact through HB 
with the molecules, both stabilization ["ordinary"] and destabilization ["surface transition"] 
behavior can occur depending on the specific parameters employed. To our knowledge this 
represents the first theoretical evidence for surface transitions for a confined binary mixture. 

A binary fluid mixture with isotropic nearest neighbor interactions is isomorphic with the 
spin 1/2 Ising model. Consider an Ising system confined between two symmetric walls which 
are D layers apart. J and Ji are the spin coupling parameters in the bulk and in the surface 
layers, respectively [p. A = '^^J'^ has a special value, A^, so that when A < A^ the phase 



mixed state is stabilized in the thin film ("ordinary" behavior) [HjIT^. In contrast, when 
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A > Ac the surface undergoes critical ordering even when the bulk is phase mixed ( "surface 



transition") [^,^. For the case of simple mixtures it has been shown that A < Ac and 
hence only ordinary behavior is predicted 0-0] • We shall show here that the introduction of 
a new interaction, the HB, could, under certain circumstances result in A > Ac, and thus 
yield complex surface behavior that is consistent with experiment. 

We consider a completely filled cubic lattice with its sites occupied by either A or 
B molecules. We model symmetric mixtures where the interactions between any two A 
molecules, and any two B molecules are identical. The molecules interact through isotropic 
nearest neighbor interactions characterized by the energy scale: x = {d/'^kBT){2eAB — ^AA — 
ess), d is the coordination number (=6), and e^j is the interaction energy between a nearest 
neighbor i — j pair. Therefore, x = 6/T* [T* = kBT/{eAB — ^aa)]- Each molecule has 
one "donor" and one "acceptor" which can participate in nearest neighbor HB. Since the 
molecules are structureless the donors and acceptors do not have prespecified locations. The 
HB interactions are described by two equilibrium constants k^A and k^B for the bonds be- 
tween an A donor and an A acceptor on different molecules, and for bonds between A and 
B particles [either one being the donor], respectively, kij = Pije~^^^^''^'^ [^. We employed 



kAA = 0.0275e^''^/^*, and k^B = 0.0134e^'^''"^*. While we have explored the phase behavior 
for a range of values of x, ^aa and kAB, we focus on this one set of parameters which yields 
closed loop phase behavior for the bulk mixture. 



In simulations of the bulk we employed periodic boundary conditions in all three 
directions. For thin films the periodic boundary conditions along the z direction were re- 
placed by two symmetric hard walls. The walls could only interact with the molecules 
through HB, i.e., ky^ = P^e~^^^''^'^ independent of the identity of the molecule. The Monte 
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Carlo simulations utilize the symmetry of the system and locate phase coexistence through 



the semigrand ensemble method with the condition that A/i = — = 0. /ij is the 
chemical potential of species i. An elementary Monte Carlo move is to change the identity 
of a randomly chosen particle. The move is accepted following the standard Metropolis 
criterion |^^. Another elementary move consists of the creation or elimination of a HB. 
The composition of the mixture [x = Na/ {Na + Nb), where Ni is the number of i particles] 
is variable and the binodal can be determined from histograms of its distribution, P{x). 
When the system is miscible P{x) has a maximum at x = 1/2. In the immiscible regime 
two maxima are observed at x = 1/2 ± Xf,, corresponding to the coexisting phases. 

In Figure |^ the binodals for the bulk and for three films of D=16 are shown. The bulk 
system displays closed loop phase behavior, a feature that is characteristic of many HB 



mixtures [|I^-|T9[. If one defines the critical temperatures as the maxima of susceptibility 
then, the upper critical solution temperature [UCST] is at T*^st — ^-94 ± 0.04, while 
the lower critical solution temperature [LCST] is at T*^^^ = 1.02 ± 0.02 [see Figure 3]. This 
result is obtained for system sizes of 8x8x8 and 16x16x16, and appears to be effectively 
independent of system size in this range. Note that a system with no HB interactions reduces 
to a standard 3-d Ising model with a critical temperature of Tq = 2.25. 

We now consider this HB mixture when it is confined between two hard, non-interacting 
walls. The mixed state is stabilized and the film binodal lies "inside" the bulk binodal. 
'^ucst = 1.80 ±0.05 and T;*^^ = 1.03 ±0.02. This behavior is in line with previous predictions 
that simple systems only show "ordinary" behavior when confined P-[^,P^. When specific 
interactions are allowed between the molecules and the walls, the phase behavior of the 
film changes qualitatively. This is the essential point of our paper. The two other curves 
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in Figure correspond to = 0.0134e^'^/^* and kw = 0.0134e^°/^*, respectively. In both 
cases the LCST is destabihzed, and the UCST is stabihzed on confinement. In the first case 
T:,,, = 1.90 ± 0.05, T:,,, = 0.98 ± 0.02, while T^,, = 1.90 ± 0.05, T^^,, = 0.99 ± 0.02 in the 
second case. 

To understand these issues better, in Figure ^ the composition in the surface layer, and in 
the center of the D = 16 films are plotted. Only composition values different from x = 1/2 
are plotted for clarity. In all cases the data from the center of the film virtually coincide 
with the bulk binodal. When the walls are neutral, the composition in the surface layer is 
closer to 0.5 than in the bulk. This is consistent with the "ordinary" transition behavior 
observed in this case. Similar behavior is observed near the UCST with the interacting 
walls. In contrast, in the vicinity of the LCST the surface is "more ordered" than the bulk, 
a signature of "surface transition" behavior. Figure ^ shows plots of susceptibility as a 
function of temperature for the middle and surface layers in the D = 16 film with k^ = 
0.0134e^^'^/^*''. The susceptibility of the middle layer tracks the bulk behavior, consistent 
with trends observed in Figure 2. In the vicinity of the LCST the surface layers show a 
distinct peak, consistent with the notion of a "surface transition". Notice also that in the 
case with strongest wall interactions, the film surface remains ordered even when the middle 
layers becomes mixed below the LCST. 

The transition from "ordinary" behavior around the UCST to the "surface transition" 
at the LCST occurs through the "extraordinary transition" , when the surface composition 
equals that in the middle of the film. This occurs at T'^ ^ 1.10 for k^, = 0.0134e(^-5/^*), and 
at T* ^ 1.60 with k^ = 0.0134e(i°-°/^*). 

A way to qualitatively understand these findings is as follows. When the walls are 
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non-interacting, the only effect of confinement is the loss of neighbors for molecules in the 
surface layer. Since this reduces the net unfavorable interaction energy of the system it 
stabilizes the single phase. In contrast, in the case of interacting walls, HB interactions 
occur between the wall and the molecules. Since each molecule has only one HB donor and 
one HB acceptor, the number of HBs between molecules within the surface layer is reduced 
due to the presence of the interacting walls. This effect is shown in Figure |^, where we 
plot the average number of molecule-molecule and molecule- wall HB contacts for the surface 
layer where ku, = 0.0134e^°''"^*. Since molecule-molecule HBs are one of the strong factors 
aiding the miscibility of these systems, and are primarily responsible for the LOST, it is 
clear that the single phase must be destabilized in the thin films. In fact, at the lowest 
temperatures in Figure ^ all donors and acceptors on molecules in the surface layer are 
occupied by the HB's with the walls. Thus, only isotropic nearest neighbor interactions 
occur between these molecules, as well as with molecules in the adjacent layer (i.e., z = 1). 
Since a bulk system with isotropic nearest neighbor interactions only possesses a UCST, to 
a zeroth approximation, the surface will remain ordered at all low temperature conditions. 

In summary, we have presented results of Monte Carlo simulations which show that the 
phase behavior of an associated fluid mixture in a thin film geometry can be qualitatively 
different from simple mixtures, which show only "ordinary" behavior. The phase behavior of 
such HB mixtures, which are defined by the balance of specific and dispersive interactions, 
can be disturbed at interfaces, leading to the occurrence of a surface transition. This finding 
strikingly rationalizes experimental results where surface transitions appear to be ubiquitous. 
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FIG. 1. Binodals for the bulk system (thick sohd hne), and for the D = 16 film with 

kw = 0.01346^*^/"^* (squares), = 0.0134e*^'^/"^* (diamonds), and without specific interactions 
with the wall, (circles). 
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FIG. 2. Compositions in the surface layer (filled symbols) and middle layers (hollow symbols) 
of the D = IQ film at the same conditions as fig. Circles - neutral walls with no interactions, 
squares and diamonds - walls with specific interactions. Squares -k^ = 0.0134e^'^/"^* , and diamonds 
-ku, = 0.0134e^-^/^* . Bulk binodal (thick sohd hne, same as fig. |l|) is shown for the reference. 
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FIG. 3. Susceptibilities X) obtained from the mean squared deviations of xa, vs. temperature 
for D = 16 film with A;^ = 0.0134e^-^/'^* . Filled diamonds - surface layers {z = 0, and z = 15), 
Hollow symbols - middle of the film {z = 7, and z = 8). Susceptibility with no surfaces and 
periodic boundary conditions is plotted for the reference (thick solid line). 
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FIG. 4. Number of hydrogen bonds per site in the D = 16 film with kw = 0.0134e^°/^* vs. T*. 
Circles- bonds between the walls and the surface layers (z = 0, and z = 15), squares-bonds within 
the surface layers, and diamonds- in-layer bonds in the middle layers {z = 7, and z = 8). 



